Quantum critical dynamics of the two-dimensional Bose gas 
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The dilute, two-dimensional Bose gas exhibits a novel regime of relaxational dynamics in the 
regime feT > \fi\ where T is the absolute temperature and fi is the chemical potential. This 
may also be interpreted as the quantum criticality of the zero density quantum critical point at 
/i = 0. We present a theory for this dynamics, to leading order in 1/ ln(A/(fcsT)), where A is 
a high energy cutoff. Although pairwise interactions between the bosons are weak at low energy 
scales, the collective dynamics are strongly coupled even when ln(A/T) is large. We argue that 
the strong-coupling effects can be isolated in an effective classical model, which is then solved 
numerically. Applications to experiments on the gap-closing transition of spin gap antiferromagnets 
in an applied field are presented. 

PACS numbers: 75.10.Jm 05.30.Jp 71.27.+a 



I 

O 

o 



> 

in 

m 
o 

in 
o 

c3 



i 

c 

o 
o 



X 



I. INTRODUCTION 

Despite the widespread recent theoretical and exper- 
imental interest in quantum phase transitions, a direct 
quantitative confrontation between theory and experi- 
ment has been difficult to achieve for systems in two 
and higher spatial dimensions. A major obstacle is that 
it is often difficult to tune system parameters over the 
range necessary to move across a quantum critical point. 
Furthermore, for many examples where such tuning is 
possible, the theory for the quantum critical point is in- 
tractable. Consequently, the analysis of the data is often 
limited to the testing of general scaling ansatzes, without 
specific quantitative theoretical predictions. 

A class of quantum phase transitions have recently 
been exceptionally well characterized in a variety of ex- 
periments. These experiments study the influence of 
a strong applied magnetic field on insulating spin-gap 
compound»i*2iMiiii&. The low lying spin excitations be- 
have like spin S z = 1 canonical Bose particles, and the 
energy required to create these bosons vanishes at a criti- 
cal field H — H c , which signals the position of a quantum 
phase transition* with dynamic critical exponent z = 2. 
In spatial dimensions d = 3, the quantum critical fluc- 
tuations are well described by the familiar Bose-Einstein 
theory of non-interacting bosons, and no sophisticated 
theory of quantum criticality is therefore necessary to in- 
terpret the experiments. The upper-critical dimension of 
the quantum critical point is d = 2, and the boson-boson 
interaction vanishes logarithmically at low momenta. So 
naively, one expects that the d — 2 case is also weakly 
coupled, and no non-trivial quantum critical behavior ob- 
tains. 

The primary purpose of this paper is to show that the 
above expectation for the quantum-criticality of d = 2 
Bose gas is incorrect. While the pairwise interactions be- 
tween the bosons are indeed weak, the collective proper- 
ties of the finite-density, thermally excited Bose gas pose 
a strong-coupling problem. We will demonstrate here 
that an effective classical model, which can be readily 



numerically simulated, provides a controlled description 
of this problem. This will allow us to present predictions 
for the evolution of the dynamic spectrum of the d = 2 
Bose gas across the quantum critical point at non-zero 
temperatures. 

Our results can be applied to two-dimensional spin- 
gap antiferromagnets, in the vicinity of the gap-closing 
transition induced by an applied magnetic. Recent 
experiments^^ on piperazinium hexachlorodicuprate 
(PHCC), (C 4 Hi 2 N 2 )Cu 2 Cl 6 , will be compared with our 
results in Section ITTT1 These experiments are able to eas- 
ily access the finite temperature quantum-critical region, 
which we will describe by our theory of the dilute Bose 
gas in d — 2; this Bose gas is in a regime of parameters 
which is typically not examined in the context of more 
conventional atomic Bose gas systems^. 

Results similar to those presented herein also apply 
to other quantum critical points with upper-critical spa- 
tial dimension d — 2 i.e. while the zero temperature 
properties can be described in a weak-coupling theory, 
all non-zero temperature observables require solution of 
a strong-coupling problem. A prominent example of such 
a quantum-critical point is the Hertz theory of the onset 
of spin-density-wave order in a Fermi liquicUi. Results 
for the thermodynamic and transport properties can be 
obtained by the methods presented herein. We maintain 
that such results are necessary for understanding exper- 
iments, and that previous theoretical resultsii*i2ii£ are 
inadequate for a quantitative analysis. 

Our results will be presented in the context of the 
phase diagram of the d — 2 Bose gas shown in Fig. fj] 
as a function of the boson chemical potential, [a, and 
temperature, T. At T = 0, there is a quantum critical 
point at [i — 0. For /i < 0, the ground state is simply the 
vacuum with no bosons, while for /i > 0, there is a finite 
density of bosons in the ground state. In the spin gap 
antiferromagnets, /i = g^B(H — H c ), where H is the ap- 
plied field, H c is the critical field, g is the gyromagnetic 
ratio, and \xb is the Bohr magneton. 

Fig^presents results in terms of contours with equiva- 
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FIG. 1: Color density plot of the value of the dimensionless 
ratio U / R as a function of the chemical potential jj, and the 
temperature T of a dilute Bose gas in two spatial dimensions. 
The energy scale U, defined in Eq. is a measure of the 
pairwise interaction between the atoms, while R, defined in 
Eq. @, is an energy scale controlling collective excitations. 
The physical properties of systems with the same value of 
U/R axe the same, apart from the change in the value of the 
energy scale R. The Kosterlitz-Thouless transition occurs at 
the contour labelled Tkt all— ^ U / R w 34, and the superfluid 
phase is present for larger U/R. The contours shown are for 
equal spacings of values of R/U, with the yellow region on the 
lower right including values R/U — > 0. Our primary results 
here are in the quantum critical region, with intermediate 
values of U/R at fi = and T > 0. 



lent physical properties, up to an overall energy scale (i?) 
whose T and /i dependence will be explicitly presented 
below, as will the equations determining the shape of 
the contours. One of the contours is the position of the 
Kosterlitz-Thouless (KT) transition of the d = 2 Bose 
gas, which is present only for /i > 0. We are primarily 
interested here in the quantum-critical region^ which is 
roughly the region with k B T > In particular, we 
will be able to quantitatively examine the signature re- 
laxation rate ~ k B T/h of the quantum-critical regimes: 
the "Bose molasses" dynamics. Our results and methods 
also describe the crossover to the KT transition, as well 
as (in principle) the region with T < Tkt- 



A. Summary of results 

We summarize here the universal aspects of our results 
of the quantum-critical Bose dynamics in two spatial di- 
mensions. The continuum quantum field theory of the 
critical point has logarithmic corrections to scaling; con- 
sequently, the properties of the continuum theory do have 



a logarithmic dependence upon a non-universal ultravio- 
let cutoff. We will show that this cutoff dependence can 
be isolated within a single parameter; all other aspects 
of the theory remain universal, and can be accurately 
computed. 

We will be interested in the continuum Bose gas theory 
with the partition function 



Z R = 



S B 



Vifi{r 1 T)e- SB/R 

h/k B T 

dr I d' z r 



7^ + |HV^| 2 -/#| 2 
ot 2m 



yH 4 



(i) 



For PHCC, the mass m can be directly determined from 
the dispersion of the S z = 1 excitation. We will present 
results primarily in the small |/i| quantum critical re- 
gion of Fig ^ although our formalism can be extended 
to other regions, including across the Kosterlitz-Thouless 
transition into the "superfluid" phase. The bare interac- 
tion between the Bose particles, Vq, is renormalized by 
repeated interactions between the particles, in the T ma- 
trix, to the value 



Vr = 



Airfr 2 



mln(A/vV + (^7 



(2) 



where A is a high energy cutoff, and the square-root func- 
tion in the argument of the logarithm is an arbitrary, con- 
venient interpolating form. The parameter A is the sole 
non-universal parameter appearing in the predictions of 
the continuum theory. For a spin gap antiferromagnet 
like PHCC, we expect A ~ J, where J is a typical ex- 
change constant. For our purposes, it is convenient to 
rescale Vr to a parameter U, which has the dimensions 
of energy: 



u ^2jnkzT VR = 



8Trk B T 



(a/\/ m 2 + (k B m 



(3) 



Our universal results for the continuum theory are 
predicated on the assumption that the logarithm in 
Eqs. J2J and is "large". At first glance, it would ap- 
pear from Eq. (|2J) that the quantum theory of the Bose 
gas is weakly coupled in the leading-logarithm approxi- 
mation. However, as will be clear from our analysis (and 
has been noted in earlier worksI^iSiiSiil) , this is not the 
case: although pairwise interactions are weak, the col- 
lective static and dynamic properties of the gas remain 
strongly coupled even when the logarithm is large. We 
shall argue, that to leading order in the logarithm, these 
strong coupling effects can be captured in an effective 
classical model. The latter model is amenable to straight- 
forward numerical simulations, and so accurate quan- 
titative predictions become possible, whose precision is 
limited only by the available computer time. Previous 



3 




20 

U/R ls 

16 
14 
12 
10 



0.05 



0.1 



0.15 —j— 1 . 0.2 

T/A 



FIG. 2: T dependence of R at the quantum critical point, 
/i = 0, obtained by solving Eq. @. A is a non- universal 
ultraviolet energy cutoff. 



FIG. 3: T dependence of the dimensionless ratio U/R at the 
quantum critical point, /i = 0. 



analyses2i2ii& of the dilute Bose gas in two dimensions 
assumed (when extended to the quantum critical region) 
that lnln(A/(fesT)) was large; a more precise version of 
this condition appears below, from which it is clear that 
this condition is essentially impossible to satisfy in prac- 
tice. We will not make any assumption on the value of 
such a double logarithm, and only make the much less 
restrictive assumption on a large single logarithm. 

The characteristic length and time scales of the 
quantum-critical Bose gas are set by a dimensionful pa- 
rameter, which we denote R. Like U, we choose this to 
have dimensions of energy; so a characteristic length is 
h/y/2mR, while a characteristic time is h/R. The value 
of R is determined by the solution of the following equa- 
tion 



R 



2tt 



(4) 



Note that this defines a R > for all — oo < fj, < oo. 

To understand the magnitude of the various scales, we 
now discuss approximate solutions of Eqs. © and Q 
at the quantum critical point, /i = 0. The estimates 
below should not be used in place of the full solutions in 
applications of our results to experiments. For the energy 
scale associated with R we obtain the estimate 



R~k B T 



41n(lln(A/(fc B r))) 
ln(A/(fcnT)) 



0. 



(5) 



Note that the energy scale R is logarithmically smaller 
than &bT: this will be the key in justifying an effective 
classical description of the dynamics. Numerically, we 
can easily go beyond Eq. JSJ , and obtain the full solution 
of Eq. (0} at jit = 0; this is shown in Fig. With the two 
dimensionful parameters, U, and R, at hand, the reader 
will not be surprised to learn that the effective classi- 
cal theory is characterized by the dimensionless coupling 
U /R. fndeed, as we will see, the ratio U/R behaves like 
an effective Ginzburg parameter for the classical theory. 
From Eqs. Q and JSJ) we estimate 



U 



2tt 



R ln(±ln(A/(/c B T))) 



/x = 0, 



(6) 



in the quantum-critical region. Again, we can go beyond 
the asymptotics, and obtain precise values for U/R by 
numerical solution of Eqs. (0J and J3J, and the result 
appears in Fig EI So, unless A/(&bT) is astronomically 
large, the ratio U/R is not small, and the classical theory 
is strongly coupled. As we have already noted, we will not 
make the assumption of large double logarithms here, but 
instead obtain numerical results for the strongly coupled 
theory. 

The primary result of this paper is that the low energy 
properties of the dilute Bose gas are universal functions 
of the ratio U / R. This key result is illustrated in Fig 
where we plot the loci of points with constant U /R, ob- 
tained from Eqs. © and (gj. The physical properties 
of systems along a fixed locus are the same, apart from 
an overall re-scaling of energy and distance scales which 
are set by the value of R. Tn most cases, the universal 
dependence on U/R can be determined by straightfor- 
ward numerical simulations. As has been shown earlier, 
in a different context, the effective classical model un- 
dergoes a Kosterlitz Thouless transition^ (see Fig^l at 
U/R « 34. The values of U/R in the quantum-critical 
region are smaller (see Eq. (JSJ) and Fig^) and the focus 
of our attention will be on these smaller values. 

With an eye to neutron scattering observations in 
PHCC, we focus in this paper on the frequency depen- 
dence of the boson Green's function, fn the antiferromag- 
net, this Green's function is proportional to the two-point 
spin correlation function in the plane orthogonal to the 
applied field i.e. the correlator of S+ and S-. The spec- 
tral density of this Green's function yields the neutron 
scattering cross-section. In particular, we will consider 
the Green's function (in imaginary time) 



h/k B T 

dr I d 2 re^" r (V>(r,T)^(0,0)); (7) 



o 



after analytic continuation to real frequency, and appli- 
cation of the fluctuation-dissipation theorem, we obtain 
the corresponding dynamic structure factor S^p(uj). One 
of our main results is that this structure factor obeys the 
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scaling form 



tuv U 
~R'' R 



(8) 



We will numerically determine the universal function 
here for a range of values of U/R in the quantum-critical 
regime. 

Our results for S^(u>) are obtained directly in real time, 
and so do not suffer ambiguities associated with analytic 
continuation. For a significant range of values of U/R 
of relevance to quantum criticality, we found that our 
numerical results could be fit quite well with the following 
simple Lorentzian functional form for xv , ( w ) 



^ ' = ~~R~ 



jR 



(hw - Rujq) 2 + (jR) 2 



(9) 



where the scaling form in Eq. JSJ implies that the dimen- 
sional numbers Z, luq, and 7 are all universal functions of 
the ratio U/R. Notice that this describes a neutron reso- 
nance at frequency Ruj /h with width R'y/h; our primary 
purpose here is to provide theoretical predictions for the 
temperature dependence of these observables. Our nu- 
merical results for the values of Z, u>o and 7 as a function 
oiU/R appear in Figs. 010 and |H1 later in the paper. 



II. QUANTUM CRITICAL THEORY 

This section will obtain the properties of the contin- 
uum theory in Eq. Q which were advertized above. 

The analysis of th e static prop erties of Eq. has 
been outlined in Ref. Il4ll5ll6ll7t The key step is the 
integrate out all the u> n 7^ modes of ip to obtain an 
effective action only for the zero frequency component. 
Among the important effects of this is to replace the bare 
interaction Vq by a renormalized interaction Vr obtained 
by summing ladder diagrams 



V 



4irh 2 



1 + (mV (47rfi 2 )) ln(A/(fc B T)) mln(A/(fc B T)) ' 

(10) 

The co-efficient of \ip\ 2 is also renormalized, as we spec- 
ify below. The resulting effective theory for the zero 
frequency component is most conveniently expressed by 
defining 



y/2mk B T 
and rescaling spatial co-ordinates by 

h 



(ii) 



(12) 



This yields the following classical partition function 



Z c = Jv9( r )e- S ' 



Sr. = 



d 2 r 



|V r *| 2 + J R|*| 2 + ^H 4 



(13) 



Here the energy U is as defined in Eq. (J3J, while the 
'mass' R is given by 



R = -fi 

2U f d 2 k 



1 



1 k 2 -fi 



.(14) 



The most important property of this expression for R is 
that the integral of k is not ultraviolet finite, and has a 
logarithmic dependence on the upper cutoff. However, 
this is not a cause for concern. The theory Z c is itself 
not a ultraviolet finite theory, and its physical properties 
do have a logarithmic dependence on the upper cutoff. 
Fortunately (indeed, as must be the case), the cutoff de- 
pendence in R above is precisely that needed to cancel 
the cutoff dependence in the correlators of Z c so that the 
final physical results are cutoff independent. This im- 
portant result is demonstrated by noting that the only 
renormalization needed to render Z c finite is a 'mass' 
renormalization from R to i?, as defined by 



R = R + 2U 



d 2 k 1 
4^2 k 2 + R' 



(15) 



Here, the R in the propagator on the r.h.s. is arbitrary, 
and is chosen for convenience. We could have chosen a 
propagator l/(fc 2 + cR) where c is an arbitrary numeri- 
cal constant; this would only redefine the meaning of the 
intermediate parameter R, but not the value of any fi- 
nal observable result. Combining Eq. I|15fl with Eq. i|14fl . 
we observe that the resulting expression for R is free of 
both ultraviolet and infrared divergences. The momen- 
tum integrals can be evaluated, and lead finally to the 
expression for R already presented in Eq. 

When expressed in terms of the renormalized param- 
eter i?, the properties of the continuum theory Z c are 
universal (i.e. independent of short-distance regulariza- 
tion). They are defined completely by the length scale 
1/VR and the dimensionless ratio U/R. The field \& is 
also dimensionless, and acquires no anomalous dimen- 
sion. This means e.g. the equal time correlations obey 
the scaling form 



s(fc) = <Mfc)| 2 ) 



1 

R 



k U 

7I'# 



(16) 



where $ is a universal function. Note that here, and in 
the remainder of the paper, we have rescaled momenta 
corresponding to Eq. (|T2*)l 



(17) 



so that k 2 has the dimensions of energy. The function $ 
can be easily computed in a perturbation theory in U / R, 
and results to order (U/R) 2 appear in Ref. 13 However, 
we are interested in values of U/R which are not small, 
and for this we have to turn to the numerical method 
described in the following subsection. 
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FIG. 4: Universal dependence of the equal-time correlation at 
zero momentum on U /R. The parameter Z is that appearing 
in Eq. ©. 



A. Numerics: equal time correlations 

We will analyze numerically Z c by placing it in a square 
lattice of spacing a, and verifying that the correlations 
measured in Monte Carlo simulations become a indepen- 
dent and universal in the limit a — ► 0. 

The partition function on the lattice is 



,-Sc 



(18) 



The parameter Rl is noi equal to the parameter i? above. 
Instead, the mapping to the quantum theory has to be 
made by requiring that the values of the renormalized R 
are the same. In the present lattice theory we have 



Rl — R 



(19) 



-2U 



71 dk x 



2tt 4 - 2 cos(k x ) - 2 cos(k y ) + Ra 2 



B. Dynamic theory 

We now extend the classical static theory above to un- 
equal time correlations by a method described in some 
detail in Refs. ITil andll^. As argued there, provided the 
scale R < ksT, the unequal time correlations can also 
be described by classical equations of motion. In the 
context of perturbation theory, the reduction to classical 
equations of motion is equivalent to the requirement that 
it is a good approximation to replace all Bose functions 
by their low energy limit: 



1 



o/T _ 1 



T 

id 



(20) 



From Eq. (JSJ we observe that the requirement on R is 
satisfied in the quantum-critical region. It also holds 
everywhere in the superfiuid phase, with fi > 0, where 
R becomes exponentially small in l/{kgT) (as can be 
shown from Eq. JI}). However, it does fail in the low T 
'spin-gap' region with fj, < 0, where R > ksT. We will 
not address this last region here, although a straightfor- 
ward perturbative analysis of the full quantum theory is 
possible here, as noted in Ref. [l9|. Some results on the 
perturbation theory appear in Appendix lAl 

The classical equations of motion obeyed by \P are 
merely the c-number representation of the Heisenberg 
equations of motion obeyed by ip. With the rescalings 
in Eqs. (fTl"1) and ljl"2|) . these are 



. fc 0* SS C 
dt <5** 



(21) 



Following the reasoning leading to Eq. (TBJ, it follows 
that the correlations of the \& evolution described by 
these equations of motion obey the scaling form 

S(k,t) = (*(fc,t)**(fc,o)) = i* t (J=, 

. (22) 

where the scaling function $ t can be determined nu- 
merically, as we describe below. After the rescalings in 
Eqs l|ll(l and 1)12(1 . we conclude that the Fourier trans- 
form of S(k, t) to S(k, uj) is related to the dynamic struc- 
ture factor S^{u)) defined below Eq. J7J by 



We can measure lengths in units of R, and for each 
value of U, determine Rl from Eq. i|19|) . and test if the 
Monte Carlo correlations are independent of a in the limit 
of small a. The resulting correlations then determine the 
scaling function $ in (|16|) . 

We used this method to determine the values of the 
function Q(0,U/R) for a sample set of values of U/R 
appropriate to the quantum-critical region. Note that 
4>(0, U/R) — Z, where Z is the amplitude appearing in 
the dynamic function in Eq. ©. We used the Wolff clus- 
ter algorithm (as described in Ref. IhH) to sample the en- 
semble specified by Z c l. Measurement of the resulting 
correlations led to the results shown in Fig|3] 



S^{u) = k B TS(0,uj) 



(23) 



We now describe our numerical computation of <i> t . We 
begin by sampling the ensemble of values specified by 
Z c l, as in the previous section. Once the ensemble is 
thermalized, we choose a typical set of values of as 
the initial condition. These are evolved forward in time 
deterministically by solving the equations of motion 



eft 



Ua 2 \*i\ 2 Vi 



(24) 



j n.n.i 



where t = t/h. For each initial condition, this defines a 
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FIG. 5: Real part of the time-dependent structure factor 
S(0, i) obtained from our numerical simulation. The results 
above are for U/R — 8.33. Two different values of R are 
shown, and the axes have been scaled to demonstrate the 
scaling collapse as required by Eq. 1221 . The fit with the 
functional form in Eq. 12611 is also shown. 
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FIG. 7: The universal dependence of the dimensionless fre- 
quency coo on U I R. The points are the results of the numeri- 
cal simulation, while the line is a best fit polynomial used for 
interpolation. 
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FIG. 6: As in Fig.0 but for the imaginary part of the corre- 
lation function. 



tfi(t). Then 



* t (p t m,u/R) = j£ 




(25) 

for a lattice of A 2 sites, and where the average is over 
the ensemble of initial conditions. Sample results from 
such simulations appear in Figs. and 

We have fit each observed time evolution to the func- 
tional form 



$ t ^Ze 



(26) 



where Z, luq, 7 are numbers determined from the fit. A 
sample fit is shown in Figs. |S] and |SJ As is clear from the 
figures, this form provides an excellent fit over a substan- 
tial time window. Taking the Fourier transform of this 
result, and using Eq. I|23|) . we obtain a result for S^(uj) 
in the form in Eq. 0, with the same parameters Z, loq, 
7- 
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FIG. 8: As in Fig Q but for the dimensionless damping 7. 
The full line at small U / R is the result of the perturbation 
theory obtained in Eq. 1A4I . 



The results for Z appeared earlier in Fig. ^ and our 
results ljq and 7 as a function of U/R appears in Figs. 
and [S] We can now combine the results in these fig- 
ures with those in Fig [31 and obtain predictions for the 
T dependence of the resonant frequeny and damping in 
the dynamic structure factor in Eq. at the quantum- 
critical point. These results appear in FiglSlandllUI 
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FIG. 9: The universal dependence of the oscillation frequency 
uioR on T at the quantum-critical point fi = 0, obtained from 
Figs. H and 
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FIG. 10: The universal dependence of the damping rate 7-R on 
T at the quantum-critical point [i — 0, obtained from Figs. 
and|U 

III. CONCLUSIONS 

This paper has argued that the quantum-critical dy- 
namics of the two-dimensional Bose gas (and of other 
quantum critical points in two spatial dimensions) rep- 
resents a strong-coupling problem. Nevertheless, an ef- 
fective classical description was obtained to leading log- 
arithmic order, allowing tractable numerical simulation. 
The primary results of this simulation at the quantum 
critical point appear in Figs. I§1 and llOl which specify the 
parameters appearing in Eq. JSJ. 

In comparing our results to experiments, we have to 
set the value of the unknown high energy cutoff A. For 
the range of parameters shown in Figs. and 1101 we 
find a roughly linear dependence of loqR and jR on 
T, with u R ps 0.11(fc B T) and jR « 0.10(fc B T). The 
experiments^*^ on the quantum-critical point of PHCC, 
also observe a roughly linear dependence of peak fre- 
quency and width on temperature, but with different co- 



efficients; the peak frequency ~ 0.62fc B T, while the width 
is pa 0.23fc B T. Choosing a different range of A for the 
theory, e.g. assuming the experiments are in the range 
T/A < 0.01, will change the theoretical predictions for 
luqR and jR, but does not improve the agreement with 
experiments. 

For the peak frequency, the work of Ref. suggests 
an origin for the above discrepancy. A similar theory 
was used in that paper to obtain the predictions of the 
peak frequency, but keeping the full lattice dispersion 
for the spin excitations, and the quantum Bose function 
values for the occupation numbers: good agreement was 
found between such a theory and the experimental ob- 
servations. The continuum and classical limits taken in 
the present paper were avoided. 

For the width of the spin excitation, we expect that 
the damping is more strongly dominated by the low en- 
ergy and low momentum excitations, and so the present 
continuum, classical theory should yield a more accu- 
rate description of the experiments. This is indeed the 
case, relative to the poor accuracy of the peak frequency 
in the continuum theory. Nevertheless, a discrepancy of 
a factor of ps 2 remains between our present quantum 
critical theory and the experimental observations. It is 
possible that taking the classical theory also on the lat- 
tice will improve agreement with experiments. However, 
there are also corrections to the present continuum the- 
ory which could improve the situation: in particular, at 
higher order in Vr, there appear renormalization of the 
time-derivative term in Eq. qTp . This renormalization 
would change the time-scale in the classical equations of 
motion, and so change the overall frequency scale of the 
results in Figs. I§l and ITUl 
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APPENDIX A: PERTUBATION THEORY 

Here we present the results of a direct perturbative 
computation of x^p- To order U 2 , perturbation theory 
yields^ 

x ^ = -fo + iLsH (A1) 

where 



E(iw) = ~2U 2 J2 J 



<fh d% 1 

4^2" 4^2" (_* ei + k{ + R)(-ie 2 + k 2 + i2)(-i(ei + e 2 ) + ^ + (h + k 2 ) 2 + R) ' 
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The frequency summation is done most easily by partial fractions and yields 

, „/ U \ 2 [d 2 k 1 d 2 k 2 [n{k 2 +R)-n((k 1 + k 2 ) 2 + R)} [n{k 2 + R) - n^h + k 2 ) 2 - k 2 )] 



\k B T J J Air 2 At: 2 ioj + (h + k 2 ) 2 - k\ - k\ - R 

Now we analytically continue to real frequencies, and take the imaginary part at w = R/h, which is the leading order 
position of the pole in Eq. i|Al|l . This yields 

TmS(R/h) = ^(J^) J ^^Hk 2 1 +R)-n((k 1 + k 2 ) 2 +R)] 

x [n(k 2 + R) - n((fei + k 2 f - k 2 )] 5{{k x + k 2 ) 2 -k\- k 2 ) 

= " 27r (^)V ^^Hkl + R)-n(k 2 + k 2 +R)] [n(k 2 ) n(k 2 + R)] 

x 5{{k! + k 2 f -k{- k\). (A2) 



Now we can do the angular integration because the angle appears only in the argument of the delta function, and 
obtain 

im£(R/n) = - y 2nkBT J J q <*iy o dk 2 ^ {k 2 +R)/[kBT) - 1)(e (fc|+i?)/( fcB T) _ _ e -(ki+ki+ R )/{k B T)) -( A3 ) 

This result is a function of R/(kgT) which has to be evaluated numerically. However, it is instructive to examine its 
value in the classical limit, upon applying the approximation in Eq. H2()(l . when we obtain 



ImE(R/h) = -—j- / dki / dk 2 
An z Jq Jo 



(k 2 + R)(k 2 + R){k 2 + k 2 + R) 



= -^y^. (A4) 

\RJ 16tt v ' 

Notice that explicit factors of k B T have dropped out, and consequently Eqs. (|A1|) and (|A4|> are consistent with the 
scaling form Eq. It is also interesting to compare the value of the classical limit in Eq. (|A4|> with that obtained 
from the full quantum expression in Eq. I)A3J| . For R = k B T/2 (which are roughly the quasiparticle energy values 
obtained in the experiments on PHCG&), Eq. l)A3j) evaluates to a value which is 16% smaller than Eq. (|A4(I . This 
gives us an estimate of the error made by representing the quantum-critical theory by classical equations of motion. 
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